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ABSTRACT 

We show that intermittency of noisclike emission, after propagation through a scattering 
medium, affects the distribution of noise in the observed correlation function. Intermittency also 
affects correlation of noise among channels of the spectrum, but leaves the average spectrum, 
average correlation function, and distribution of noise among channels of the spectrum unchanged. 
Pulsars are examples of such sources: intermittent and affected by interstellar propagation. We 
assume that the source emits Gaussian white noise, modulated by a time-envelope. Propagation 
convolves the resulting time series with an impulse-response function that represents effects of 
dispersion, scattering, and absorption. We assume that this propagation kernel is shorter than the 
time for an observer to accumulate a single spectrum. We show that rapidly-varying intermittent 
emission tends to concentrate noise near the central lag of the correlation function. We derive 
mathematical expressions for this effect, in terms of the time envelope and the propagation 
kernel. We present examples, discuss effects of background noise, and compare our results with 
observations. 

Subject headings: methods: data analysis - techniques 



1. INTRODUCTION 

The observed spectrum of a source with intrinsically intermittent emission, including propagation effects, 
is the focus of this paper. Propagation effects can often be described as a convolution, in time, with an 
impulse-response function or "propagation kernel". This tends to smear out any intrinsic variation of the 
source. In particular, if the emission from the source is noiselike, but modulated in time, then propagation 
will mix emission at different times, so that the observer perceives noiselike emission of nearly constant 
intensity. However, we find that imprints of intermittency remain, under certain simple assumptions. 

Our work is motivated by observations of pulsars at meter and decimeter wavelengths. Pulsars show 
a rich variety of time variability, on timescales as short as nanoseconds, as we discuss below. Propagation 
broadens these rapid intrinsic variations with propagation kernels, with typical time spans of microseconds 
to milliseconds. The forms of these kernels vary over longer timescales, of seconds or more. Our goal is to 
distinguish between intrinsic and propagation effects, as simply and with as few approximations as possible. 

A vast variety of propagation effects can be described as a convolution of the emission, in time, with 
a propagation kernel. Such a convolution introduces correlations among samples, which an observer can 
measure. We suppose that the source emits Gaussian-distributed noise, with time-varying amplitude. Most 
astrophysical sources emit noise of this type, and it provides the standard model for pulsar emission. Because 
our assumptions are relatively broadly based, we anticipate that our results may have wider applicability. An 
observer measures spectra or correlation functions from the time series, and estimates an average, regarded 
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as the determinstic part, or signal; and variations about that average, usually regarded as noise. We suppose 
that the observer accumulates data for a single spectrum over a timescale short compared with that for the 
propagation kernel to vary, and makes statistical comparisons among spectra. 

Using our model approximations, we calculate the averages and noise for spectra and correlation func- 
tions. Among the various factors that contribute to noise, we focus in particular on self-noise, produced by 
the noiselike emission of the source. We find that intermittency introduces correlations of the noise among 
spectral channels in the spectrum, and alters the distribution of noise in the correlation function. The aver- 
age spectrum and correlation function of the source are the same as those of a non-varying source, with the 
same propagation kernel. Intermittency of the source affects only noise, but affects it in a way that can be 
calculated. Thus, noise carries information about the source that is lacking in averages. 



1.1. Outline 



In this paper, we discuss the effect of intermittent emission of a noiselike source on the spectrum and 
correlation function, in the presence of propagation effects. In the following Sj2j we introduce treatment of 
propagation as a convolution, and quantify the approximations used for this treatment. We discuss noiselike 
emission and relate it to the noise and averages of the observables. We describe observations of the pulsar 
B0834+06 as an example, and evaluate our approximations. 

In Sj3j we present formal calculations of the averages, or signal; and of noise, for the spectrum and 
correlation function. This section forms the core of this paper. We introduce notation and fundamentals of 



noise, and introduce the time-envelope and the propagation kernel, and the observed time series (£3.1 1. We 
transform to the frequency domain of the power spectrum, and present expressions for the average spectrum, 



and for the variances and covariances of noise of spectral channels (£3.2). We then transform back to the 



time-lag domain of the correlation function and present expressions for the average correlation function and 



noise in the correlation function ([3.4). All of these quantities depend on the propagation kernel. However, 



only noise in the correlation function, and correlation of noise between spectral channels, depend on the 
time-envelope of emission at the source. We present examples of spectrum and correlation function for 
simple cases of source intermittency, and simple propagation kernels, in j|3.5| 



In f|4] we discuss comparisons and extensions of our mathematical theory. We present connections to 
the amplitude- modulated noise theory of Rickett ( 1975 1, describe effects of background sky or receiver noise, 



and compare our results with observations. In q5]wc summarize our results. 



THEORETICAL BACKGROUND 



2.1. Propagation as a Convolution 



A convolution in time describes effects of wave propagation through a medium with spatially-varying 
index of refraction, static in time. Multiplication of the spectrum by a frequency-dependent complex scalar 



provides an equivalent description, via the convolution theorem for Fourier transforms (Bracewell 2000) 



Applications include cell-phone communication, seismology, ocean sound propagation, radar and microwave 
circuits, adaptive optics, and interstellar radio- wave propagation ( Simon et al.||200f| Moustakis et al.pOOO 



|Bullen fc Bolt|1985|[Ziolnek|1985|[Shearer|2009[pontgomery, Dicke fc Purcell|1987HGoodman|2 000; GwETet] 



al. 1998 ). For electromagnetic waves, this fact stems from the superposition principle. A localized assemblage 
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of charges and currents with 4- vector current density J^(t,x), gives rise to a 4- vector potential A^{t,x). If 
current density varies at a single frequency ui, then is a solution to the Hclmholtz equation (iLandau & 



|Lifshitz|1975 1. In this equation, a spatially- varying wavenumber describes a time- independent, but spatially 
non- uniform, index of refraction (Ersoy 20071. The solution for A^^^x) for an arbitrary arrangement of 



refracting material exists and is unique under quite general assumptions. Helmholtz equations of this form 



arise for propagation of many kinds of waves, and a host of techniques provides means of solving them (Flatte 
1979| |Ziomekfl985{ |Goodmanf2000| |Ersoy||2007[ and references therein). 



For a particular frequency, source configuration, and arrangement of refracting material, the electric 
field at the observer equals that at a fiducial point near the source, times a complex propagation factor g. 
The propagation factor can depend on frequency. For time-independent refraction, g(u>) is independent of 
the phase or amplitude of the electric field at the source. More generally, a matrix g(uS) relates different 
source polarization states, and different source locations, to the electric field in different observed polarization 
states and observer locations. This is the S-matrix or scattering matrix, originally developed to describe 



quantum-mechanical scattering 


( Cushing 


19901, and then ; 


(Montgomery, Dicke & Purcell 


1987| Moustakis et al.||2000 



leads to a superposition of received electric fields; under some circumstances, such as a sufficiently small 
source, the details of the source distribution are unimportant ( Gwinn et al.||1998 Goodman 2000; Simon et 
aLl[200T ). 

Multiplication by the scalar g(uj) in the frequency domain becomes convolution by its Fourier transform 



g{n) in the time domain, where k is the time lag (Bracewell 2000 see also eq. 13 below). Indeed, for any 
linear time-invariant system, of arbitrary dimensionality and complication, one time series input and one 



time series output are always related by convolution with such a kernel (Rabiner & Gold||1975 Oppenheim, 



Willsky, & Nawab 1997). The kernel g(n) is known as the "propagation kernel" or the "impulse-response 



function" . Its width is the typical scattering time r, expressed in terms of time lag k. 



2.2. Spectrum and Dynamic Spectrum 

2.2.1. One spectrum 

This paper is concerned with statistics of spectra and correlation functions. In this paper, we use the 
terms "one power spectrum" , "one realization of the power spectrum" , or "an estimate of the power spec- 
trum" for the square modulus of the discrete Fourier transform of a time interval of data, of N v samples 
of duration 5t each. The Fourier transform of the power spectrum is the autocorrelation function. In- 
terferometry combines two time series: the product of the Fourier transform of one with the conjugated 
Fourier transform of the other is one realization of the cross-power spectrum. The Fourier transform of the 
cross-power spectrum is the cross-correlation function ( Thompson, Moran, fc Swenson|2001 Perley, Schwab, 
& Bridle 19891. In this paper, we use the term "spectrum" to refer to either the power spectrum or the 



cross-power spectrum. Sometimes the power spectrum is known as the "autocorrelation spectrum" or sim- 
ply the "auto spectrum" ; and the cross-power spectrum is sometimes known as the "cross spectrum" . The 



Fourier transform has both discrete and integral forms; Thompson, Moran, & Swenson (20011 describe the 



two forms and their relations. This paper is concerned with finite spans of sampled data, and hence we use 
the discrete Fourier transform. In some formal contexts, one realization of the power spectrum is known 
as the "periodogram" , whereas the "spectrum" is the integral Fourier transform of a continuous time series 
( Broersen||2006 ) ; however, we follow traditional usage here. 
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We call an average of spectra over an ensemble of statistically-identical spectra with different realizations 
of noise "the ensemble-average spectrum" or simply the "average spectrum" . We denote specifically averages 
over time, frequency, or propagation kernel, when we invoke them. Effects of intermittency or propagation 
are not included in our statistical ensemble: only those of noise. Averaging the autocorrelation function over 
an ensemble of realizations of noise yields the "average autocorrelation function" . 

This paper deals with scattering that is approximately time-independent, over the time interval for 
accumulation of a single spectrum, f the medium remains time-independent for many accumulation times, 
consecutive spectra are statistically identical, and the observer can average them together to increase signal- 
to-noise ratio. We discuss these considerations in detail in f|3] 



2.2.2. Dynamic spectrum 



Over longer times, the propagation kernel varies and the spectra differ. This paper is not directly 
concerned with this regime, although we briefly discuss analysis of data with such variations in Q A simple 
way of dealing with time series that vary with both frequency and time is the dynamic spectrum or, more 
formally, short-time Fourier transform (STFT) ( |Bracewell||2000| |Chen k Ling||2002"l |Boashash||2003[ ) . The 
observer divides a longer time series into segments of N v elements each, and forms one spectrum from 
each segment. Individual spectra can be averaged together for periods shorter than the timescale for the 
propagation kernel to change. The time series of such spectra is the dynamic spectrum, a two-dimensional 
function of time and frequency. 

The 2D Fourier transform of the dynamic spectrum is the autocorrelation function in the lag-rate 
domain. The Fourier transform of the spectrum is the autocorrelation function, as noted above. The Fourier 
transform along the time axis converts the time series of autocorrelation functions to frequency, known as 
"rate" . The square modulus of the resulting lag-rate autocorrelation function is the "secondary spectrum" . 
The secondary spectrum is particularly well-suited to study of the largest delays produced by interstellar 
scattering, which is observed to take the form of "scintillation arcs" in this domain (Hill et al. 2003). 
Theoretical description of scattering is quite advanced, and makes extensive use of the dynamic spectrum 
and secondary spectrum to infer the statistics of scattering material, and distribution in particular cases 
flRickett||l977l |Walker et aI1[2004l [COTdes et al.|[2006|. 



2.3. Approximations 

In this section we relate our fundamental approximations to the properties of source, scattering material, 
and instrumentation for observations. The assumption of time-independent refraction is an approximation 
for all the applications mentioned above; but in many cases the approximation is excellent. These are 
fundamentally instrumental approximations in that their accuracy depends on the time period the observer 
samples. If the time for accumulation of a single realization of the spectrum is short, the medium can be 
treated as time- independent. We also require that the time for accumulation of a spectrum is longer than 
the time span of the propagation kernel. This is likewise an instrumental constraint. Finally, we assume 
that the source is noiselike: more precisely, the source emits a random electric field, drawn from a Gaussian 
distribution at each instant and uncorrelated in time. The variance of the Gaussian distribution gives the 
amplitude of the source. This paper deals specifically with the effects of intermittency on short timescales, 
where the variance of the distribution changes during the time for accumulation of one spectrum. 
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2.3.1. Time variation of propagation kernel 



Propagation can be treated as time-independent when the time to accumulate a spectrum is shorter than 
the characteristic time for the propagation medium to change significantly, at the retarded time. Propagation 
includes dispersion and scattering. Dispersion changes only slowly with observing epoch ( Backer et al.|1993 l, 
far too slowly for variations to affect a single spectrum. Scattering changes much more rapidly. For most 



media, the propagation kernel for scattering has a relatively rapid rise, with a much slower fall (Williamson 



1972 Boldyrev fc Gwinn|2003 ) . Within that envelope, phase and amplitude vary rapidly as the electric field 
along different paths cancels or reinforces. The propagation kernel changes with time, with changes in the 
scattering medium or locations of source and observer. 

Motion of the material in the propagation medium produces quasistatic evolution of the solution to the 
Hclmholtz equation. A slow variation of the propagation factor g with time represents the effect ( Ziomek 



1985). For small-angle forward scattering by angle 6 and speed v, the local rate of change of phase is about 
duv/c, corresponding to a time for change by about f radian of order t& k, c/{6vv). For observations of 
galactic objects at distances up to a few kiloparsecs, we expect scattering material to move at v ~ 10 to 
100 km sec -1 , and 9 some fraction of the observed scattering angle. In many cases, motion of the source sets 
the timescale for scintillation; the scattering material is assumed "frozen" as it moves across the line of sight 
Goodman 2000). This assumption is the foundation of velocity measurement of pulsars, via scintillation 



flGupta|1995| ). In this case, the expression for td is similar, with v now the speed of the source perpendicular 



to the line of sight, and with factors of order 1. This assumption is also fundamental to most theories of 
material responsible for scattering (see Walker et al.|2004 1 . For galactic sources observed at meter or shorter 
wavelengths, typical timescales td are seconds or longer. The typical accumulation time for an individual 
spectrum N„St is about a millisecond. Thus, the approximation is quite good. We present a specific example 



in i 2.5 below 



2.3.2. Doppler shift 



The approximation of time-independent refraction holds when a single spectrum has spectral resolution 
insufficient to detect Doppler shifts from motions of the source or medium. In this case, a propagation 
factor g(ui) describes scattering accurately. For small-angle forward scattering, the change in frequency from 
Doppler shift is 8v ~ Qvvjc. This expression is the same as that in the previous section, and indeed the net 
effect of Doppler shifts can be expressed as a time derivative of g(oj) (Ziomek 1985). For typical galactic 
10 -11 for scintillation arcs, and much less for the typical scattering (Walker et al. 2004 1. 



sources. 



Av/v 



The frequency resolution of a typical spectrum, 1/N v 8t, is typically a few kilohertz. Thus, these Doppler 
shifts are not detectable in individual spectra, at observing frequencies less than 100 GHz. Effects of the 
Doppler shifts do appear in dynamic spectra, as quasi-static evolution of the propagation kernel in time. We 
present an example in |2.5| 



2.3.3. Time span of the propagation kernel 

In our calculations, we require that the time span of the propagation kernel, r, is shorter than the time 
for accumulation of a single spectrum, N v 8t. The kernel includes effects of both scattering and dispersion. 
For scattering time and scintillation bandwidth Ai^ = l/2wTd, this becomes the requirement that the 
spectral resolution be finer than Avd- This requirement is mildly violated for some observations of scattering 
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at the largest r, as we discuss in S2.5 



In practice, the timescale for dispersion tjjm can be longer than the timescale for the spectrum, as we 
discuss in §2.5| This does not present a fundamental obstacle, because the contribution of dispersion to the 
propagation kernel can be inverted. Dispersion correction via multiplication by the conjugate of the expected 
phase of g(u>), is precise but computationally intensive and is usually reserved for the most interesting cases 
(Hankins 1971). Incoherent dedispersion alleviates effects of dispersion (Voute et al. 2002), but is not a 
pure deconvolution and can be expected to leave some artifacts. Ideal measurements would involve coherent 
dedispersion, or long accumulation times for single spectra, or both. 



2.4. Noiselike Emission 



We assume that the emission from the source is noiselike. Specifically, we assume that the source emits 
Gaussian white noise: the electric field at each instant is random and drawn from a Gaussian distribution, and 
the field at different instants is uncorrelated ( Papoulis]|l991 ). Many sources, including astrophysical sources 
under almost all conditions, emit such radiation, because they comprise many superposed, independent 
radiators (Dicke 1946 Evans et al.||1972" ). The distribution of electric field is then Gaussian because of the 
central limit theorem. We do not assume that the emission of the source is stationary; in other words, we 
suppose that the variance of the underlying Gaussian distribution may vary with time. 



The noise in an observable V is the departure of a given measurement from its statistical average: 

SV = V-(V) n . 



(1) 



This average is calculated over an ensemble of statistically-identical realizations of self-noise and background 
noise, denoted by the subscripted angular brackets (...)„. "Self-noise" or "source noise" is the contribution 
of the noiselike emission of the source to noise of observables ( Kulkarn"i||1989 Anantharamaiah et al.fl991 



Gwinn|2004 2006). We quantify this effect in £3.3.2 below, and describe observations in 



Gwinn et al. 



(2011) 



Background noise includes Gaussian noise from sky, ground and antenna emission, and noise in the receiving 
system. The ensemble does not include variations in the flux density of the source or of the propagation 
kernel. We treat these effects with the functions / and g, defined below. Consequently, in this paper, we do 
not approximate the statistical ensemble by averages over time or frequency, and we do not include effects of 
changes in time envelope or propagation kernel on periods longer than the period to form a single spectrum. 
The variance of the observable V quantifies that noise: 



(5V 2 ) n = 6V 2 ee (V 2 ) n (V)l 



(2) 



If the distribution of noise is Gaussian, then V is completely characterized by its mean, the signal; and its 
variance, the noise in V (|Papoulisp9l"]). 



2.5. Example Observations 

In this section we describe observations, and evaluate the approximations discussed in §2.3[ We consider 



observations of scintillations of the pulsar B0834+06. Gwinn et al. (2011) report single-dish observations at 
Arecibo and very-long baseline interferometry observations on the Arecibo-Jodrell Bank baseline; |Brisken| 



et al. (2010) report very-long baseline interferometry observations using a network of antennas, including 
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the Arecibo-Green Bank baseline. Both observed at frequency v 
Av f 



320 MHz, with observing bandwidth of 



MHz. Both used incoherent dedispersion (Voute et al. 2002). Instrumental parameters important 



for the approximations are as follows: The time for accumulation of a single spectrum was N v 8t — 0.33 msec 
for single-dish and 0.128 msec for VLBI observations of Gwinn et al.; and 8 msec for Brisken et al. All of 
the observations averaged spectra together in time to increase signal-to-noise ratio: Gwinn et al. averaged 
over Aig = 10 sec, and Brisken et al. over 5 pulses, or At$ = 6.4 sec. All observations gated synchronously 
with the pulsar pulse; thus, only a few individual spectra of the many in these time intervals contributed to 
averages. 

Parameters of the source were measured in these and previous observations. Dispersion broadens a 
short pulse by tdm ~ 23 msec at v = 320 MHz (Lorimer & Kramer 2004). Scintillation produces strong 
variations of intensity across the spectrum, characterized by a typical scintillation bandwidth Avd = 0.57 
MHz ( Gwinn et al.|20lT |. The characteristic timescale of the propagation kernel is = 1/2ttAv ( i — 0.3 /xsec. 



The time span of the kernel is a few times greater. The scintillation timescale is = 290 sec (Gwinn et 



al. 20111. This is the timescale for variation of the propagation kernel. Note that these observational and 



instrumental parameters meet with the requirements in { 2.3 duration of the propagation kernel is less than 
accumulation time for a spectrum, which is less than the averaging time, which is less than time for variation 
of the propagation kernel: < N u 8t < At$ < t^. Note that dispersion violates the requirement that the 
propagation kernel is shorter than the accumulation time for a spectrum: tdm > N v St. However, incoherent 
dedispersion mitigates the violation. The inferred angular broadening is about 9d ~ 1 mas, so scattering is 
small-angle forward scattering. 

Fine-scale structure appears in the dynamic spectrum of pulsar B0834+06, corresponding to a scintilla- 
tion arc ( Hill et al.|2003 1. The arc represents about 3% of the spectral power, so it is a small, but interesting, 
effect of propagation ( Gwinn et al.||2011 |. The arc extends to delay as great as r a = 1.2 msec, and to rates 
as high as us a = 2n x 50 mHz (Brisken et al. 2010). This rate corresponds to a timescale for variation of 
the propagation kernel of t a ~ 20 sec . For the arcs, the approximations hold well for the observations of 
Brisken et al.: < N v 8t < At s < t a . Gwinn et al. miss the longest-timescale part of the propagation kernel 
for the arcs, but this contains relatively little spectral power. The inferred angular deflections range up to 
6 a i=s 30 mas, still small-angle forward scattering. 



3. MATHEMATICAL THEORY 

In this section, we calculate average spectrum and correlation function, and noise, using our approxi- 
mations of |2.3| Table [T] summarizes the symbols we use to describe emission, propagation, and observation. 
Table [2] summarizes our mathematical results. 



3.1. Definitions and Notation 

3.1.1. Noise 

We model source emission as Gaussian white noise of constant amplitude. This noise is multiplied 
by a time-envelope function to model intermittent emission. This noiselike, but nonstationary emission is 
then subjected to spectral variation, which expresses the effects of propagation. An instrument converts an 
electric field to a complex time series via a Hilbert transform 



BraceweiT||2000 1 and then converts a limited 



- 8 - 



band of frequencies to a baseband time series (Lorimer & Kramer 2004). The observer forms products of 



samples of the baseband time-series to form spectra or correlation functions. 

The spectra may vary greatly with time; as an example, intrinsic emission from pulsars varies greatly 
with time, and interstellar propagation imposes time-varying spectral modulation, by convolution with the 
propagation kernel. It is essential that the observer not normalize the spectra, or otherwise impose a time- 
varying gain. This would change the variances of the underlying contributions to noise, independently of the 
variations introduced by source emission and the particular realization of propagation, so that underlying 
contributions would be impossible to disentangle. Similarly, we define normalizations of theoretical quantities 
in this section. 

Noiselikc emission has an electric field ej. Here, the index j represents time. Each sample ej is drawn 
from a complex Gaussian distribution with zero mean and unit variance. We suppose that the underlying 
noiselike emission is "white": different samples are uncorrelated ( Papoulis||1991 ). Then, 

( e J e l)n = S jk . (3) 

Here, Sjk is the Kronecker delta symbol, which is 1 if j = k and otherwise. Because the noiselike emission 
has no intrinsic phase, 

(e i e fc )„ = (e*e* k ) n = 0, (4) 

whether j = k or not ( Gwinn|2006 ). We leave open the possibility of different noiselike sources: for example 
source and background noise, or background noise for two different antennas; we denote these different 



electric fields with superscripts: e\ 



and so on. Different "noises" are uncorrelated: 



Sjk S x 



(5) 



We discuss the consequences for background noise in [4.2 



3.1.2. Flux Density 

The intensity of the source (or, equivalently, the flux density) is proportional to the mean square of 
the electric field. For a noiselike electric field, intensity is the variance 
constant flux density has stationary electric field. Variations in the time envelope and in the propagation 
kernel can both change the flux density. 

To describe sources of arbitrary flux density, we add a prefactor of the average flux density, I s . Because 
we consider only pairwise combinations of sources of noise in this mathematical section, we can omit the 
prefactor I s in this mathematical section without loss of generality. We normalize the noise using Eq. [5] 
above, the time envelope using Eq. |6j and the propagation kernel using Eq. |J effects of any flux density 
variations are subsumed into I s . We reintroduce I a in Qfor comparison with other work and observations, 
where multiple components of noise with different variances must be taken into account. This parameter 
also takes into account the possibility of variations in average flux density among spectra because of changes 
in the time-envelope or propagation kernel. 



Dicke||1946 ) . A noiselike source of 



3.1.3. Time Envelope 



The emission may be intermittent: the electric field may be modulated in time. In this case, the electric 
field is a white, Gaussian random variable; but it is not stationary (Papoulis 19911. We accommodate this 
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by scaling ej by the emission envelope fj . This function expresses a change in the standard deviation of the 
emission, without changing its noiselike character. Multiplication by fj in the time domain is equivalent to 
convolution with its Fourier transform in the spectral domain; we make much use of this fact below. We 
demand that the time-envelope of emission, fj, not change the average intensity: 

E uu = (6) 

3=0 

Without this normalization, variations in source emission will change the means and variances of the cor- 



relation function and spectrum, making impossible the calculation of signal and noise ([3.3 and [3.4). We 
assume that fj is real and positive: a change in phase of fj has no effect on the statistics of the underlying 
noiselike emission. Noise from different sources can have different time-envelopes f x , f v : noise from a source 
may vary, while noise from a receiver remains constant. For later convenience, we introduce the products of 
time-envelopes: 

# = % = fjfi P7 = ftf v i (7) 

If the time envelopes are the same, /J = /J, all the /3's are the same: (3* = 0% = (3 xy . 



3.1.4- Propagation Kernel 
After emission, propagation affects the intermittent, noiselike emission by convolution in the time do- 



main, as discussed in [2.1 above. We parametrize this effect by convolving with the kernel g^. This kernel 
includes all propagation effects: dispersion, scattering, and even spectrally-varying absorption. It can also 
include the spectral response of the observer's instrument. Convolution with gk is equivalent to multiplica- 
tion by its Fourier transform in the spectral domain. Similarly to fj, we demand that the propagation kernel 
not affect average intensity: 

E I- (8) 

3=0 

Like Eq. [6j this normalization allows calculation of signal and noise in §3.3| and j ]3.4| As argued below, 
products of <?'s are the average correlation function or spectrum, a or p. For a stationary white-noise source, 
without time or spectral variations or propagation effects, fj — 1, and g X j — <5jo- 



3.1.5. Time Domain and Wrap Assumption 

To make a measurement, an observer forms products of two, possibly distinct, electric fields xi and yi . 
To obtain these fields, the underlying noiselike emission from the source e is multiplied by / and convolved 
with g: 

N u -x 

"*= E Sift-tf-i W 

3=0 

and similarly for y. Here, the span of the observation is N v samples. The mean intensity is (xx*) n , and 
mean interferometric visibility is (xy*) n . More generally, the time series can be correlated with different 
lags and Fourier transformed to form a spectrum; or they can be Fourier transformed and then multiplied 
together to form a spectrum ( Thompson, Moran, fc Swenson||2001 ). 
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For convenience, and for compatibility with the circularity of the discrete Fourier transform in £ 11 we 
make the "wrap" assumption throughout this paper ( Gwinn|2006 1 : 



X, 



— Xi, 



(10) 



for any quantity X. This simplifies calculations, and accurately describes an "FX" correlator (Chikada et 



al7||1984 1 . Some "XF" correlators violate this assumption, but effects on noise are calculable and relatively 
minor (Gwinn et al. 2011b). Moreover, we form circular convolutions and correlations (Rabiner & Gold 



1975 Oppenheim, Willsky, & Nawab 1997). Padding of the N v samples with zeros recovers the common 



case of ignoring the wrap. 

The different fields xi and yi might arise at different observing stations as in interferometry, or from 
different sources (such as source noise and background noise, or different spatial regions of one source). 
For interferometry of a single pointlike source, the underlying noiselike emission e and time-envelopes / 
are identical for two stations x and y, but propagation g may be different. We will retain the identifying 
superscripts in the results throughout this section, so as to extend the model to include multiple components, 



for discussion in S4.2 



3.2. Frequency Domain 

3.2.1. Fourier transform, Parseval's theorem, convolution theorem 

The Fourier transform relates the electric field in the time domain Xg, with that in the frequency domain 
Xk- We adopt the normalization convention for the Fourier transform and the inverse transform: 



Xh 



iV„/2-l 

e X t 

t=-N v /2 

iV„/2-l 



(11) 



Xi = — E ^ ke X k 



N v 



k=-N u /2 



Here, N v is the number of spectral channels. This is the discrete form of the Fourier transform (Thompson, 
Moran, fc SwensonpOOlj ) . 



Parseval's Theorem states that power is conserved by the Fourier transform. Thus, 

N„-l JV„-1 



E W = n v e 



(12) 



fc=0 



Power is conserved separately for signal and noise; thus, this equation holds for the square of statistically- 
averaged (V) n and for noise power SV 2 as well as for individual time series, spectra and correlation functions. 

The convolution theorem, sometimes known as the Faltung theorem, states that the Fourier transform 
of the convolution of two functions is the product of their Fourier transforms. Thus, 



V k U k = E e * 

rn— 



N v -1 



E V « U > 



m—£ 



(13) 



1=0 
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for any functions V and U. Using the fact that the Fourier transform of X^ e is X£ (see Eq. 11 ) 

Nv-l 



N u —X 

ij^km i/*n . 

j=0 



(14) 



m— 

For the inverse transform, the convolution theorem takes the form: 



k-j 



(15) 



fc=0 



3=0 



3.2.2. Data series in the spectral domain 

In the spectral domain, we deal with the Fourier transforms of the time series. These Fourier transforms 
are the spectral series Xk, fjk- These spectral series, in turn, can be described in terms of the Fourier 
transforms of the underlying noise, e|, e\\ the time-envelopes /£, f£; and the propagation kernels g\. 
Using the definition of xi and the convolution theorem Eq. |13| we find 



Xk 



N v -1 
3=0 



Because the Fourier transform of white noise is white noise (Papoulis 1991), {e^e v t ) 



>h-e 



(16) 



6 e . Here, we 



introduce the notation S e to express a value of 1 if e x and e y are identical, and if they are different. Because 
of Parseval's Theorem, J2k fkfk* = N v , an d J2k 9k9k* = 1> an< ^ similarly for the quantities depending on y. 
Because the time-envelope functions are real, = f x * k . 



3.3. Spectrum 

3.3.1. Average Spectrum 

When two spectral series are multiplied together, XkHk, their product is one realization of the cross- 



power spectrum, as in an FX correlator (Chikada et al. 1984 Thompson, Moran, & Swenson 2001). This 
is 



1 



N„ 



rk 



XkVk 



zv* 

-m m ' 



j^r E St fk-t g ? 9k fk- 

V e,m=Q 

If a single series is multiplied by itself, its square modulus is one realization of the power spectrum: 

1 



®k = ~^~ x kXk, 



(17) 



(18) 



The expression for a k in terms of g x , f x , and e x is identical to Eq 

I i < ■ .ui.iLi.iiiin Kii I 



17 



but with y — > x. The expression for 



Of. is analogous, but with x — > y. The factor of 1/N V is required for self-consistent normalization. 



The statistically-averaged spectrum pk (or ctk for the power spectrum) is the deterministic part of this 
realization. We calculate this as an average over noise, while holding the time-envelope and propagation 
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kernel constant. We again denote the average with subscripted angular brackets 

N v -1 



p k = (r k ) n = ~g%~9t*^- E fk-Jf-m (efe»T)„. (19) 

v i,m=0 

where we have made use of the fact that different samples of noise are uncorrelated, Eq. [5] Again, the 
symbol S e indicates a factor of 1 if the noise series e x and e y are identical, and if they are different. This 
expression defines the ensemble-average cross-power spectrum pk- Interestingly, this average depends only 
on the propagation kernel g. If the underlying noise e is different for x and y (as for combination of source 
noise with background noise) then pk = 0, as the symbol 8 e indicates. The autocorrelation functions also 
depends only on the propagation kernel: 



~9l~9l\ (20) 



and likewise for a k . Because of the normalization of the propagation factor g discussed in £3.2.2 the 
normalization of a is X)fc=o 1 ^ k = Because of Parseval's Theorem, for our normalization of the Fourier 
transform, the flux density of the source averaged over time is equal to the flux density summed over the 



spectrum. In our calculations we have set both to 1 (see £3.1.2). Although these results are presented in 



the context of spectral variations imposed by scattering, they are valid for a noiselike source with any sort 



of imposed spectral variations, including line absorption. The notation here is identical to that of Gwinn 



3.3.2. Noise in the spectrum 

Noise in the spectrum is the departure of a realization of the spectrum from the average. The variance 
in a spectral channel characterizes noise in that channel. More generally, we can calculate the covariance of 
noise in two (possibly different) channels as 

S (hr2) = {fkffln ~ {fk)n{f*t)n (21) 
1 N„-l 

- J_ V 7, x f x h x * f x * /p x p x *\ h y f v n v * f v * /pVpV*\ 

- jy2 9kJk-p9e li-r \ e p e r ) n 9lH~s9k Jk-q \ e s e q /„ 

^ N v -1 N v -1 

= a k a T vT X! fk-p fp-l JT- X! $T-t fk-q 
V p=0 V q=0 

- n x n y * R x R y * 

- a k a e Pk-lPk-t 

We have used the fact that K = —flk> as weu as * ne ever-crucial fact that samples of the emitted electric 
field are uncorrelated, Eqs. 

H and U 

Here, /3k is the convolution of the time-envelope function with itself, in 

the spectral domain: 

v p=0 

and similarly for /3 k . The convolution theorem shows that f3 is the Fourier transform of f3, Eq. [7] Because of 
Parseval's Theorem, and our normalization of / in the time domain, /3q = M = 1. This expression remains 
unchanged if the underlying noiselike time series e k and e v k are uncorrelated: then the mean correlation of x 
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and y is pk = as noted in { 3.3.1 above, but the noise is the same, given by Eq. 21 For the power spectrum 
of xe, the same expression gives the noise, with of course ffc — > af. on the left-hand side of the equation, 
and with a k and all for the same station and propagation kernel on the right: y — > x. The analogous 
construction yields noise for a v k . 



For a single spectral channel k, Eq. 21 simplifies greatly, and becomes: 



'kin 



a k al 



(23) 



This equation implies that intermittent emission, as expressed by the time-envelope function, does not change 
the noise in an individual spectral channel at all. In fact, this equation states that the signal-to-noise ratio 
in a single spectral channel is 1: 8{rkf k ) /\pk\ 2 = 1 (see Broersen 2006). This is consistent with the usual 



radiometer equation (Dicke 1946), if one observes that the number of samples per spectral channel is 1 . If the 



observer averages together N t independent samples of the spectrum, then the signal-to-noise ratio increases, 



as the noise falls proportionately to 1/N t , or more slowly if amplitude variations are present (Gwinn et al 



2011) 



In general, the time-envelope function can affect the covariance of noise between spectral channels, as 
Eq. [21] suggests. Indeed, the normalized covariance of noise between two channels depends only on the 
time-envelope function and the separation of channels (fc — £): 



J k-e 



PI* 



(24) 



The normalized covariance does not depend on the propagation kernel or the particular channels k or I. For 
a source emitting white noise of constant intensity, j3k = 5k , an d the covariance of noise in different channels 
is zero. For a source emitting a single emission spike during the integration, /3k = 1, and the normalized 
covariance of noise is 1 for all pairs of channels: the noise is perfectly correlated. We discuss these examples 
further in 93.51 below. 



The power spectrum is real, so the noise must real as well; but the cross-power spectrum can be complex. 
In this case, the noise forms an elliptical Gaussian distribution in the complex plane, and characterization 
of the noise requires another variance. This measures the degree of elongation of the distribution in the 
complex plane: 



5(f k fe) 



1 

N?. 



AT„-1 

E 

p,q,r,s— 



Pi) 



9kfk- P ~9Tf!-s {K^)jl}t-r9 V k*f V k- q (^ef> r 



(25) 



~ nxy 

Pk Pi Pk_ t 



B xy * ■ S 



Here, we have made use of the hybrid quantity /3^ v , the Fourier transform of (3 xy in Eq. [7] Eqs . 



25 



and 



21 



show that noise is strongest in phase with signal (see |Gwinn 2006 Eq. 20; Gwinn et al. 2011| §4). If the 
underlying noiselikc series ej: and e v k are uncorrelated, S(fkfi) = 0; in that case, the distribution of noise in 
the complex plane is a circular Gaussian distribution. 
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3.4. Lag Domain: Correlation Function 

3.4-1- Average correlation function 

Correlation characterizes the observed time series in the lag domain, Fourier conjugate to the spectral 
domain. The cross-correlation function r K is given by: 

1 Nu-l 
v e=o 

We define the autocorrelation function a% analogously, but of course by auto-correlation of identical time 
series. We can determine the statistical averages through a direct average of Eq. [26] or by using the 
convolution theorem for the Fourier transforms of xg and yt with Eq. |19| We find 

N u -1 

P K = {r K )n = E St&n-te. (27) 
1=0 
AT„-1 

< = (On = E 9i9tX K , 
1=0 

and similarly for a y K . Like the average spectrum, the average correlation function depends only on the 
propagation kernel, not on the time-envelope function. Note that our normalization condition (Eq. [8| 
implies ckq — 1. Again, if the underlying noiselike time series ef and e\ are different, then p K = 0. 



3.4.2. Noise in the correlation function 



Now we consider the noise in the lag domain. We can obtain this by the inverse Fourier transform of 



the noise in the spectral domain, Eq. 21 



<5( r «<) 



K 1 Kin V K/n\i K /n 
1 N "~ 1 

k,£=0 
Nu-1 



(28) 



N?, 



-y 



N u -1 



x 

or ol 



Thus, noise in the lag domain is the correlation of the time-envelopes at the two stations, convolved with 
the product of the autocorrelation functions at the two stations. Here, we have used Eq. [14] to re-express 
ftk-l Pk—ii an d the facts that /3 is real, and that a K = a*_ K . The expression for the noise in the autocorrelation 
function, 5(a K a*), is identical, with y — > x. As this equation shows, in the lag domain, the noise in a single 
lag is affected by the time-envelope: unlike signal in either domain, or noise in one channel of the spectral 
domain. If the underlying noiselike series ef, ef are distinct for x and y, then the noise is still given by Eq. 
|28| unchanged. For the power spectrum of the time series xg, the same expression again gives the noise, with 
of course r K — ► a% on the left-had side of the equation, and with a K and (3 K all for station x on the right: 
y — > x. The analogous expressions hold for yg. 

The practitioner may be surprised that intermittent emission affects noise in the lag domain, but not in 
the spectral domain. After all, noise in one domain is the Fourier transform of noise in the other. However, 
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we calculate the variance of noise; and correlations in one domain affect variances in the other, a consequence 
of the convolution theorem. The order of imposition of time variation and convolution is important: if the 
emission from the source were first convolved with a propagation kernel, and then modulated in time, then 
the time-envelope would affect noise in the spectral domain. 

In general, correlation is complex in the lag domain: p K and a K are complex. We require a second 
variance to fully characterize it. This can be found by Fourier transform of Eq. |25| 

S(r K r K ) = (r K r K ) n - (r K ) n (r K )n (29) 
= 4i E e-^e-^ p k ft fa ft* ■ S e 

v k,t=0 

= E e- i ^ k+ ^~p k p e e- i ^ k -^^y^ +v -S e 

v k,e,fi,iy=0 

JV„-1 iV„-l 



This gives the ellipticity of the distribution of noise in the complex plane. Here, we have again used Eq. |14| 
If the underlying noiselike series ef, e\ are distinct for x and y, then 5(r K r K ) = 0, and the distribution of 
noise is circular. The expression for the noise in the autocorrelation function of the time series X£ is again 
identical, with r K — > a x K and p->aas well as y — > x. 



3.5. Examples 

In this section we present some simple examples of the theory above. In the examples, we suppose that 
xi and yi are identical, so that our observations measure intensity for a single noise source. We therefore 
consider autocorrelation functions a K and spectrum a k , and their statistically-averaged counterparts a K and 



3.5.1. Constant-intensity source 

In the simple and traditional case of constant source intensity over the integration time, source electric 
field is simply Gaussian white noise, with constant variance. Such emission characterizes most continuum 
sources. Propagation introduces correlations between samples at different times, without changing the 
variance of the distribution. Figure [T] shows such a convolution, for a short span of random electric field and 
a schematic propagation kernel. The electric field is drawn from a complex Gaussian distribution at each 
sample. Real and imaginary parts have variances of so that the intensity is 1. The propagation kernel is 
a one-sided exponential, gj — tjt exp(—t/T s ) for t > 0, with characteristic time T s — 5. This kernel satisfies 
the normalization condition for gj, Eq. [8] This resembles effects of scattering, with a rapid rise and slow fall, 
although the complicated amplitude and phase variations are absent. As the figure shows, the convolution 
smooths the electric field. Because the kernel is normalized, the intensity is unchanged. From this time 
series, one can calculate a single realization of the autocorrelation function. Alternatively, one can Fourier 
transform the electric field and form its square modulus: this is one realization of the spectrum. 

The average autocorrelation function and average spectrum agree with predictions of |3.3.1| and §3.4.1| 
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Figure [2] compares results of numerical simulation with the theoretical results. For a source with white 
noiselike emission with constant intensity, the time-envelope is constant: fg = 1, and so fig — 1 and /3& = 5ko- 
The average autocorrelation function a K is simply the autocorrelation function of the propagation kernel gj , 
as Eq. 27 indicates. In this case, it is a two-sided exponential, a K = exp(— \k\/T s ). The average spectrum is 
the square modulus of the Fourier transform of Qj, as Eq. 19 indicates: a Lorentzian with half- width at half 
maximum N u /(2ttT s ), for our discrete Fourier transform. In the examples in Figure [T] and [2j T s approaches 
N u , so the forms depart slightly from these analytic forms in our calculations and simulations. 

In our model, noise arises purely from source noise. Noise agrees well with theoretical estimates, as 
Figure [2] shows. For a constant-intensity source, the noise in one channel of the spectrum is equal to the 
square modulus of the average spectrum in that channel (Eq. 23 Brocrscn 2006 ) . The signal-to-noise ratio 
is 1, but this can be increased by averaging together samples of the spectrum with the same propagation 
kernel. The noise of the autocorrelation function is constant as a function of lag. Eq. [28] shows that it is 
S(a K a* K ) — X)^o 1 \ a "\ 2 ~ 2/T s - For the individual realizations, noise must be identical at positive and 
negative lags: this is a consequence of the fact that the power spectrum is always real, so that a K = a*_ K , 
even for an single realization. 



3.5.2. Single- spike emission 

As an example from the opposite extreme, suppose that the time envelope consists of a single spike: a 
delta- function. The time-envelope is fj = ^N v Sjq. This satisfies Eq. [6j the normalization condition for fj. 
The electric field at the source is nonzero only for one sample, with the random electric field of that sample 
drawn from a Gaussian distribution. The observed time series is the convolution of the resulting spike with 
the propagation kernel, yielding a copy of the propagation kernel gj, scaled by the random electric field of 
the spike, and shifted in time. The observed time series is a single, one-sided exponential, with random 
amplitude and phase. Figure [I] shows the resulting time series, and effects of the propagation kernel. 

For our propagation kernel, the autocorrelation function is a two-sided exponential, times a random 
gain factor |eo| 2 ; and the spectrum its Fourier transform. However, averaging over many such realizations 
recovers the same averages as in the constant-intensity example, as Figure [T] shows. 

Noise is quite different. The noise is a single gain-like factor, which changes the amplitude of the 
spectrum but does not change its form. Noise in the spectrum turns out to be simply the square of the 
average spectrum, just as in the case of white noise. However, each realization of the spectrum is a copy 
of the average spectrum, scaled by a random factor, so that noise is perfectly correlated among channels. 
Similarly, in the lag domain, each realization of the autocorrelation function is a copy of the average, scaled 
by some amount. Consequently, the variance of noise in the lag domain is the square of the autocorrelation 
function. The distribution of noise is not uniform, as for a constant-intensity source. Noise is concentrated 
near zero lag, and falls rapidly away from it. Figure [2] shows this example as well. 



3.5.3. Multiple spikes 

Emission in multiple, randomly spaced spikes is intermediate between continuous emission and a single 
spike. Cordes ( 1976 ) proposed such shot noise as a fundamental model for pulsar emission. We suppose that 
the pulsar emits n e single-sample spikes during the observing period of N v samples, and that emission is 
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zero at other times. We suppose that all spikes have equal amplitude, and that they are randomly spaced. 
Then f p = ^/ N v /n e or 0. We consider a propagation kernel with phase as well as amplitude variations. 

Figure [3] shows an example. The source emits 6 spikes of equal average intensity. The time-envelope fg 
is shown at the bottom of the top panel of the figure. Each spike has a random electric field, drawn from a 
Gaussian distribution. The spikes are then dispersed by a propagation kernel gf. a complex sinusoid with 
a declining envelope. We suppose that x and y are identical: we measure a K and We average results 
over many realizations with different noise, but with identical time-envelope and propagation kernel. The 
averages are computed with spikes at identical locations (that is, fg is fixed), but varying source noise eg. 
As the upper panel shows, convolution with the kernel g disperses the spikes in time, and allows them to 
interfere if they overlap. The form of the kernel is apparent. The average correlation function in the middle 
panel shows the imprint of phase variations in the propagation kernel, and of its characteristic decline with 
time. Any effect of correlation between different spikes is absent in the average, because differences in phase 
of eg between spikes cancel it. The correlation function displays the required normalization ao — 1 and 
symmetry a K = a*_ K . 

Correlation between different spikes contributes to noise, as shown in the lowest panel. Noise in the 
correlation function is largest at the central lag, but nonzero to the highest lags. Vertical lines at the bottom 
of the lower panel show the autocorrelation of the time envelope, /3k- The distribution of noise displays the 
overall behavior of j3 K . The spiky j3 K is smoothed by convolution with the square of the average correlation 
function a K . This leads to a relatively smooth "floor" level of noise for all k. This "floor" arises from Eq. 
[28} and is the same for real and imaginary parts. Eq. [29] dictates dependence of the noise on phase, and is 
nonzero only where individual convolved spikes interact with themselves via the terms p K + M and p K -u- This 
takes place at k — ¥ and k — > N v /2. Consequently, the noise depends on phase in these regions. The "bump" 
at large positive and negative lag arises from the latter. If desired, it can be removed by zero-padding the 
time series. 



3.5-4- Random emission events 

In general, we expect spikes, or other patterns of emission, to vary in location among different samples 
of the spectrum. They may also vary in number and strength, and have finite width and various shapes. 
Description of these effects requires averaging over the time-envelope fg. These effects can be included in our 
model by dividing a much longer span of observations into segments, and finding effects of time integration 
by convolution of the spectrum. Here, we briefly discuss the form of such spectra, and effects of averaging 
over time-envelope as well as noise. 

Because the average correlation function a K is affected only by the propagation kernel gg, variations of 
the time-envelope fg do not affect it. They affect only the noise, through the autocorrelation function of j3 K 
in the second sum in Eq. |28| Suppose that the emission consists of a number n e of randomly-placed spikes 
of width one sample, with uniform variance, as above. Then (3 K — N u /n e or 0. The correlation function of 
/3 K , times 1/N%, is a peak of area N% at k — 0, and n e (n e — 1) peaks of area N^/n^ distributed elsewhere. 
The central peak results from correlation of each spike with itself, and the other peaks from correlation of 
spikes at different times. Convolution with the square of a K smooths the spiky correlation function of /3 K to 
a more nearly flat "floor" level, as Figure [3] illustrates. The area of the central peak is the product of areas 
of the convolved functions, A = J2 K l a ft| 2 ) an d area of each of the outlying peaks is Ajn\. 

Averaging over many spectra with the same number of spikes at different locations, and the same 
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propagation kernel, leaves the central peak unchanged, while smoothing the noise floor to a constant level. 
The central peak contains noise power A, and the floor contains total power A(n e — \)/n e , or power per lag 
of A(n e — l)/(n e N v ). For a single emission event n e — > 1, the peak dominates and the noise floor vanishes, 
as expected (see j |3.5.2 ). For many spikes n e » N„, the noise floor dominates and the noise is distributed 



approximately uniformly in lag, as expected for stationary noiselike emission (see ^3.5.1). As Cordes (1976) 
notes, superposition of sufficiently many spikes results in stationary noise. 

If the emission takes the form of "events" of some typical width w e , then the resulting average auto- 
correlation function takes the same form: a central, broadened peak and a surrounding noise floor. The 
normalization of /3 K demands that the central peak have unit area (although now width ~ w e ), and that the 
outlying peaks again each have area N v /n e (and width ~ w e ). Noise in the autocorrelation function takes a 
similar form: a central peak of noise power A, and a uniform noise floor of power per lag A(n e — l)/(n e N u ). 

Figure [4] shows an example of random emission events: the average autocorrelation function and noise 
for a time series containing 8 randomly-placed time- envelopes of length 13 and constant height, and an 
exponentially-declining real propagation kernel. The average includes random placement of the events, as 
well as varying noise. Signal reflects only the propagation kernel; noise reflects propagation kernel and the 
time-envelope. The noise shows a floor of the expected level. The central noise peak is broadened by the 
autocorrelation function of the individual emission events. In contrast to results for 1 or more spikes in 
Figures [2] and [3j the noise peak is broader than the signal peak. This is because of the greater width of the 
emission events w e , and consequent broadening of (3 K . 



DISCUSSION 



4.1. Connection with Amplitude-Modulated Noise Theory 



Our model for pulsar emission is based on the amplitude-modulated noise theory of Rickett ( 1975 ). This 



work provides a test of whether pulsar emission is stationary Gaussian noise, multiplied by a variable time 
envelope (Rickett 1975). This model has become a standard benchmark for studies of statistics of pulsar 



emission ( |Jenet, Anderson, fc Pri"nce||2001| [Smits et al.|[2003) |Popov et al.||2009[ |Jessner et aL||2010| |. Both 
our and Rickett 's analyses are concerned with the second and fourth moments of the electric field. As a 
test of this property, iRickett (1975) invokes the correlation function of intensity Ri(k) — (I(t)I(t + K)}t, n - 



Here, the angular brackets include both an average over time, and over realizations of noise, as indicated 
by the subscripts. This function is related to the statistics described in the present paper, particularly to 
noise, the fourth moment of electric field. Rickett includes in his analysis a convolution in time, similar 
to our propagation kernel, to describe effects of a bandpass filter in the receiver. However, he does not 



include propagation effects explicitly. Smits et al. ( 2003 ) point out that propagation effects can be difficult 



to remove, and residual ones can strongly affect the observed form of Rj. 



Our treatment differs from Rickett (1975) in that we include propagation effects explicitly, including 



those with timescale t longer than the timescale of variation of the source. Rickett assumes that the minimum 
timescale for variation of the time-envelope (his w, our w e for width of an emission event) must be much 
longer than timescale associated with the kernel (his 1/B for the inverse of the receiver bandwidth, our r for 
the propagation kernel). For pulsar B0834+06, for example, this demands that flux density variations of the 
pulsar be longer than 1.2 msec, much longer than many of the variations that are observed. Our expressions 
in f|3] hold for arbitrary time-envelope and propagation kernel; of course, our kernel can include receiver 
response as well. We also include both auto- and cross-correlation, and we distinguish between signal and 
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noise contributions. We do not invoke an average over time, as Rickett does, but only over the statistics of 
the noiselike electric field of the source. 

Rickett's paper and the present work have fundamentally different aims: Rickett proposes a test of 
whether a source emits amplitude-modulated Gaussian noise in the absence of propagation effects, whereas we 
calculate effects of propagation under the assumption of amplitude-modulated noise. It would be interesting 
to test whether the source emission is noiselike, in the presence of propagation effects. This requires additional 
knowledge of, or assumptions about, propagation and source emission. In a promising approach, |Jessner| 



et al. (2010) compare the form of i?j for different pulse phases; specifically, for giant pulses from the Crab 
pulsar. Interstellar propagation is not likely to produce differences between pulse phases. Comparison of 
signal and noise for the autocorrelation functions, for different pulse phases for this object, would produce 
observables that we could compare with the present work. 



4.2. Background Noise 



Background noise often dominates self-noise and is always present. Background noise is not coherent 
with source noise. Thus, we can take the electric field as a superposition of time series x\ for the source and 
a;™ for noise, and calculate the signal and noise for each pair of fields separately. To maintain generality, 
we must relax the assumption that the time series have unit variance. We assign variances / s and /„; these 
are the average flux densities of of signal and noise, respectively. We maintain the convention that a and 



p are normalized as above, and adopt I s and I n as prefactors, as discussed in j 3.1.2 For interferometry, 
background noise is uncorrelated between antennas, and we must distinguish between backgrounds at the 
two antennas I n A and I n B- Usually background noise is of constant intensity, so that <5£ = 1 for all fc; 
generalization to spectrally-varying background noise is straightforward. In the spectral domain, we recover 
the expressions for interferometric visibility V k : 



04) „ 

s(v k v k *) 

5(V k V k ) 



(30) 



(I s a s k ) + (I s a s k ) {(I nA + I nB )} + {J„ A /„ B } 



For a relatively short baseline, p k ~ a k , and the expressions for variance become a second-order polynomial 
in the average visibilty (Vk} n for the real part, and a first-order polynomial for the imaginary part (see Gwinn 
et al. |201l| Eq. 8). Similarly, for single-dish observations, we obtain in the spectral domain the intensity I k : 



(h)n = I s a s k 



In 



(31) 



{half 



(I s a s k ){2I n } + {ll} 



Again, the expression for noise is a second-order polynomial in the average intensity (see Gwinn et al.|2011 



Eq. 11). The polynomial is a perfect square, for intensity. When averaged over N t samples in time, with 
constant flux density in each spectral channel, each term declines by a factor of 1/JVj; the polynomial remains 
a perfect square. 

If the average flux density of the source changes slowly as spectra are averaged, the polynomial coeffi- 
cients can decline more slowly than 1/N t . We present a detailed model in Gwinn et al. (2011), §2.2.2. In 



the simplest cases, such as a pulsar with a rectangular pulse, the effect can be treated as a reduction in the 
number of independent samples. For more complicated amplitude variations, the coefficients depend on the 
moments of the distribution of mean intensity, and the polynomial is no longer a perfect square. 
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4.3. Example Observations 



We propose that variations of / on timescales shorter than the accumulation of a single spectrum explain 
the apparent lack of noise at large lags for pulsar B0834+06 ( Gwinn et al.|20lT ). The distribution of noise in 
the spectrum is in agreement with that for a constant-intensity source, following the polynomial forms given 
by Eqs. 30 and 31 with corrections for amplitude variations ( Gwinn et al.|20lT §2.2.2). This noise, must be 
conserved in a transform to the correlation function, in accord with Parseval's Theorem. In this domain, the 
noise is uniformly distributed for interferometric observations, which are dominated by background noise. 
However, noise appears deficient at large lags of the correlation function for single-dish observations, which 
are dominated by self-noise. We suggest that this apparent deficiency is precisely the concentration of noise 
at small lags predicted for an intermittent source, as seen in Figures [3] and [3] and discussed in §3.5.4| 

We can form a quantitative estimate of the degree of source intermittency using the results of §3.5.4| In 
the single-dish observations, we measure a noise level of about 40% of that expected for stationary noiselike 
emission, from noise in the spectral domain and Parseval's Theorem, at lags greater than 40 /isec. About a 
quarter of this is sky or instrumental noise, which is uniformly distributed. Thus, we suggest that self-noise 
is reduced to about 30% of its constant-intensity value at large lags, by intermittency. 



In terms of the model from \ 3.5.4[ reduction to 30% indicates about 1.5 emission events per integration 
period of 300 /isec. The inferred duration of the events is w e < 40 /isec. Peaks of emission of this short 
duration, repeating at about this rate, are indeed observed for this pulsar ( Kardashev et al.|1978 ). Of course, 
alternative models, for example involving occasional powerful narrow emission events interspersed with more 
typical constant emission, can explain the observations equally well. 



SUMMARY 



In this paper, we consider the effects of intermittent emission at the source on observations of individual 
realizations of the spectrum and correlation function. We argue in §2.1| that effects of propagation can be 
described as a convolution in time, or multiplication in frequency, to high accuracy. We assume that the 
underlying emission at the source is white noise drawn from a stationary Gaussian distribution, modulated 
by a time-envelope function f(t), as discussed in [2.4 and f|3j We suppose that this emission is convolved 
with a propagation kernel g, and that the resulting electric field is correlated at the observer, as described 
mathematically in j |3.1.5[ 

We exploit the properties of Gaussian noise to find that variations of the time-envelope do not affect the 
average spectrum, the average correlation function, or the noise in individual channels in the spectrum, in 
j}3~3land [3XT] These are all determined only by the propagation kernel and the flux density of the source. 
The time-envelope does affect noise in the correlation function, and the correlation of noise between different 
channels in the spectrum. We present expressions for the distribution of signal and noise, and spectral 
correlations of noise. Table [2] summarizes our mathematical results. 

The distribution of noise in the correlation function depends on the form of the intermittency. For an 



source of constant intensity, the noise is constant with lag ([ 3.5.1 1. For a single spike, the emission is zero 



outside a central peak ([3.5.2). For multiple randomly-occurring emission events, the noise has a constant 
"floor" surrounding a central peak ( §13.5.4 1. The width of the peak is set by the timescale of the propagation 
kernel and the widths of the events. We contrast our calculation with the amplitude-modulated noise theory 



of Rickett (1975), which does not include propagation effects. 
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Background noise adds a constant offset, and a linear term, to the polynomial describing noise as a 
function of source intensity (\ 4.2 1. We apply our theory to recent observations of noise in the secondary 
spectrum of the scintillating pulsar B0834+06, and show that these suggest that the emission of the pulsar is 
intermittent on timescales shorter than 40 /isec 4.3 1 . Our theory suggests, interestingly, that the existence 
of intermittency can be inferred even when temporal smearing, via convolution with the propagation kernel, 
has reduced the degree of modulation. Additional work, now in progress, suggests that more extensive 
analysis can recover both the propagation kernel and the form of the emission, at least in statistical form. 



We thank the National Science Foundation (AST-1008865) for financial support. 
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Fig. 1. — Effects of convolution with a propagation kernel (top panel) on a time series consisting of 
constant-intensity noiselike emission (middle pair of panels), and on a single spike (lower pair of panels). 
Constant-intensity emission is drawn from a Gaussian distribution with unit average power. The spike is 
drawn from the same distribution, but with amplitude adjusted so that the mean square power is the same 
in the two cases. Circles show real parts, crosses imaginary parts. 
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Autocorrelation Function Spectrum 




-10 10 -10 10 

Lag t Channel k 

Fig. 2. — Correlation functions and spectra of constant-intensity noiselike emission and single-spike emission, 
convolved with a one-sided exponential. Solid lines show expected forms, and circles (for constant intensity) 
and crosses (for spike) show signal and variance of noise for 10 5 simulations like that in Figure [l] Upper 
panels: signals in correlation (left) and spectral (right) domains. Lower panels: noise. Left panel shows 
the correlation function, with flat curve for constant-intensity noiselike emission and sharp peak for a delta- 
function spike. 




Fig. 3. — Autocorrelation function and noise for random spiky emission and more complicated propagation 
kernel. Top: Typical time series. Vertical lines at lower edge show 6 spikes of the time-envelope f p , 
displaced downward. Solid line shows real part, and dotted line shows imaginary part, of one time series 
X£, including effects of multiplication by random eg and convolution with propagation kernel g p . Middle: 
Average correlation function for time series, with f p and g p as in top panel. Lower: Noise for correlation 
function. Solid line shows variance of real part, dotted line shows variance of imaginary part. Vertical lines 
at lower edge show /3 K , displaced downward. 
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Fig. 4. — Autocorrelation function and noise for wider bursts of source emission, with varying times of 
bursts. Top: Typical example of time series. Horizontal lines at lower edge show time-envelope of emission 
events, where f p > 0. During each of 8 bursts, width w = 13, source emission is noiselike with constant 
variance, and is convolved with the one-sided exponential propagation kernel. Solid line shows real part, 
and dotted line shows imaginary part, of one typical time series xg. Middle: Average correlation function, 
for average over source noise and over times of emission events. Lower: Variance of noise for correlation 
function. The offset from zero is the noise "floor" described in the text. Noise peak is broadened, relative 
to the peak in the middle panel, by the finite width of the bursts. 
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Tabic 1. Dictionary of Symbols 



Symbol 



Definition 



Reference Notes 



Time Series 

Xf_ 

e l 

£X 

h 
at 



Spectral Domain 



8{a k a k ) 



Pk 

s (?kfl), S(f k f k ) 

I 

V 

Is 
In 

Lag Domain 

a K 
n J 

S(a K a* K ), S(a K a K ) 
r K 

Pk 



Observed Electric Field 
Underlying Stationary Noise 
Time-Envelope Function 
Propagation Kernel 
Square Modulus of ft 

Power Spectrum 
Average Power Spectrum 
Noise in Power Spectrum 
Cross-power Spectrum 
Average Cross-Power Spectrum 
Noise in Cross-Power Spectrum 
Intensity 

Interferometric Visibility 
Average Flux Density of Source 
Flux Density of Noise 



Autocorrelation Function 
Average Autocorrelation Function 
Noise in Autocorrelation Function 
Cross-correlation Function 
Average Cross-correlation Function 
Noise in Cross-correlation Function 



Eq. | 



S 3.1.1 



3.1.3 



S3. 1.4 



Eq~j?[ 



Eq. 
Eq. 



18 



20 



3.3.2 



17 



Eq. 
Eq. 
EqsTpll 25 



19 



3.1.2 



4.2 



(4.2 



3.4.1 



Eq. 



2G 



3.4.2 



2(i 



Eq. 
Eq. 
Eqs. 281 29 



27 



a,b 
a,b 
a,b 
a,b 
a,b 



a For quantities for companion field ye, substitute x — > y. 
b Fourier transforms of these quantities indicated by accent 
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Table 2. Summary of Equations 



Description 



Equation" 



Reference 



Time-Envelope Function 
Normalization of / 
Normalization of / 
Definition of j3 
Normalization of (3 
Normalization of j3 

Propagation Kernel 
Normalization of g 
Normalization of g 

Spectrum 
Average Cross-Power Spectrum 6 
Average Power Spectrum 
Normalization of a 
Noise in Single ChanneP 
Covariance of Noise in 2 Channels' 3 
Conjugate b ' d 

Correlation Function 
Average Cross-Correlation Function 6 
Average Autocorrelation Function 
Normalization of a K 
Noise in a Single Lag c 
Conjugate 6, c 



PI = ft ft 

k = ~pi = i 

T,t 9e 9e* = 1 

Pk = 9k 9t* ■ $e 

at = 9k 9k* 
Efc«fe = 1 
S(r k if) = a%a v k * _ 

*(r*#?) = a%af_PZ_ 
S(f k ft) = PkPiPl-t.Pt 



• Pk-l 
X y* _ g 



P* = T,t 9e 9 V t *+K ■ § e 
a% = J2i9f9e+ K 
a = 1 

S(r*r* K ) = ^ E M E„ PI Pl+u 

S(r K r K ) = Jf7 E M PK+fi Pn-fi Tliv Pv V P^+v 



Eq. 



i 3.2.2 



Eq. M 
Eq. M 



^3.3.2 



Eq. | 



^3.2.2 




see ^3.2.2 
Eq. E| 
Eq.|2T" 



Eq. 
Eq. 



f 3.4.1 



Eq. 
Eq. 



a All sums run from to N v — 1. 

b If noise series e x , e y are identical, then S e = 1; if noise series are different, then S e = 0. 

c For expressions for noise in power spectrum or autocorrelation function a%, substitute y — > x, 
p — > a. 

d Vanishes for autocorrelation function a and power spectrum a. 



